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Abstract 



' The combined effective low energy QCD Lagrangians of Nambu - Jona- 

Lasinio (NJL) and 't Hooft are supplemented with eight-quark interactions. 



This work is a follow-up of recent findings, namely (i) the six quark flavour 



O 

determinant 't Hooft term destabilizes the NJL vacuum, (ii) the addition of 
a chiral invariant eight-fermion contact term renders the ground state of the 
theory globally stable; (iii) stability constrains the values of coupling con- 
stants of the model, meaning that even in the presence of eight-quark forces 
the system can be unstable in a certain parameter region. In the present 
work we study a phenomenological output of eight-quark interactions con- 
sidering the mass spectra of pseudoscalar and scalar mesons. Mixing angles 
are obtained and their equivalence to the two angle approach is derived. We 
^ ■ show that the masses of pseudoscalars are almost neutral to the eight-quark 

forces. The only marked effect of the second order in the SU(3) breaking is 
found in the r] — 7/ system. The scalars are more sensitive to the eight-quark 
interactions. A strong repulsion between the singlet-octet members is the 
reason for the obtained low mass of the a state within the model considered. 
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1. Introduction 



The fundamental fields of QCD, quarks and gluons, are unobservable dy- 
namical variables. Instead, at low energies, one observes hadrons. The most 
direct way to study their properties is the method of effective Lagrangians 
written in terms of the matter fields describing mesons or baryons. Such the- 
ories, under some circumstances, can be developed to the advanced level of 
an effective field theory. A well-known example is chiral perturbation theory 
PP, where the Lagrangian of light pseudoscalar mesons is both a derivative 
and light current quark mass expansion around the asymmetric ground state 
which is assumed to be stable. This stability is a phenomenological fact; the 
underlying theory must explain it. 

It is commonly accepted that the eight approximate Goldstone bosons 
7T, K and rj are a signal for the spontaneous breakdown of the chiral symmetry 
of QCD, which is realized in the ideal world of massless u, d, s, quarks. It is 
not excluded that effective four-quark interactions of the Nambu and Jona- 
Lasinio (NJL) type [2j are responsible for the formation of a stable chiral 
asymmetric vacuum giving a crude insight into the structure of the ground 
state of QCD 0. 

One might ask if higher order quark interactions are of importance. For 
instance, on lines suggested by an instanton-gas model, it can be argued jl] 
that there exists an infinite set of multi-quark terms in the effective quark La- 
grangian starting from the NJL four-quark interactions. The famous 't Hooft 
determinantal interaction [3] automatically appears if one keeps only the zero 
mode contribution in the mode expansion of the effective Lagrangian. This 
2Nf multi-quark term (Nf being the number of quark flavours) manifestly 
violates the Ua(1) axial symmetry of the QCD Lagrangian, offering a way 
out of the Ua(1) problem. 

Let us recall the case with the lightest flavour singlet pseudoscalar rj', 
which was for a while a deep theoretical problem, known as the £7a(1) puzzle. 
The general solution [S| has shown that the rj', being a quark-antiquark state, 
is strongly connected to the gluon world and that the C/a(1) axial anomaly is 
the reason for the rj — ir,K,ri splitting observed in nature. These conclusions 
are based only on the Ward identities and the 1/N C expansion of QCD (where 
N c is the number of colours). 

The same question has been also studied in the framework of an effec- 
tive Lagrangian which includes the meson fields and the topological charge 
density Q(x) jjj. After the elimination of the field Q(x) by means of its 
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classical equation of motion, one obtains an effective mesonic Lagrangian. It 
has been shown by Rosenzweig, Schechter and Trahern jH] that the 't Hooft 
type determinantal interaction, written in terms of mesonic fields, appears 
as the first term in the expansion which results from eliminating Q(x). From 
the Lagrangian of the model one realizes again that there are no valid the- 
oretical objections against the idea that the 't Hooft interaction and higher 
order multi-quark terms are actually present in the QCD vacuum. Part of 
these interactions have been utilized in jHj. 

Thus, it is tempting to consider the intuitive picture that describes the 
QCD vacuum with basis on a series of multi-quark interactions reflecting 
several tractable features of QCD, which include aspects of chiral symmetry 
and of the 1/N C expansion. The bosonization of quark degrees of freedom 
leads then to the desirable effective Lagrangian with matter fields and a 
stable chiral asymmetric vacuum. 

This idea is not new. The NJL type model with the Ua{V) axial symme- 
try breaking by the 't Hooft determinant (in the following we will use the 
abbreviation N JLH for this model) has been studied in the mean field approx- 
imation ^OJ- CHI for a long time. Numerous phenomenological applications 
show that the results of such an approach meet expectations. Nevertheless 
in this picture, there is an apparent problem: the mean field potential is 
unbounded from below, and the 't Hooft term is the direct source of such 
an instability (see, for instance, Eq. (3.16) in ^S])- A consistent approach 
requires obviously a stable hadronic vacuum in which the pions would live 
forever in the ideal world with only strong interactions. 

The functional integral bosonization of the model exposes new shortcom- 
ings: the system of stationary phase equations used to estimate the generat- 
ing functional of the theory Z, has several real solutions ^7j which contribute 
independently, i.e., Z = Zi + Z 2 + . . - + Z n , where n is the total number of such 
real solutions. Since only one of them (let us assume Z±, for definiteness) 
leads, at leading order, to the mean field potential, Vmf, the semiclassical 
potential V, corresponding to Z, differs from Vmf- It has been shown in J7| 
that V is also unbounded from below. Thus, we must accept that the NJLH 
model suffers from a ground state problem. 

Recently it has been argued [TH] that eight-quark interactions, added to 
the NJLH Lagrangian, might resolve the problem. Indeed, the mean field 
potential of the modified theory, Vmf, has a globally stable minimum. The 
just mentioned controversy concerning the results obtained by the mean field 
method and the functional integral approach is also removed: one can prove 
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that V = Vmf-i i-e., the number of admissible real solutions to the stationary 
phase equations can be constrained to one due to eight-quark terms. 

There is a natural question. If the eight-quark forces are so important 
for the formation of the ground state, what are the other phenomenological 
consequences of such interactions? 

In this paper we consider the main characteristics of light pseudoscalar 
mesons (J PC = h ): their masses and weak decay constants. After that 
we switch to scalars, calculating masses of the J PC = ++ quark- ant iquark 
nonet. The structure of scalars is a subject of many studies nowadays. The 
question is so complicated that it would be too naive to think that eight- 
quark forces are a panacea for the mass spectrum problem. Our aim is only 
to demonstrate the tendency. Once we understand what is changed by the 
new interactions considered in the description of the meson properties within 
the model, we can clarify the role of eight-quark forces for low-energy QCD. 

To study this matter one should choose an appropriate approximation. 
The bosonization of six- and eight-quark interactions cannot be done ex- 
actly. We will use the stationary phase method to replace the multi-quark 
vertices by purely mesonic ones and by Yukawa type interactions of quarks 
with mesons. This is a standard approach jTT], [IH]-|2I]- The subsequent 
integration over quarks is a straightforward calculation, because one deals 
here with a Gaussian integral. To obtain the effective mesonic Lagrangian 
and extract masses, we shall expand the real part of the quark determinant 
in a heat kernel series [221 123! • The techniques which are particularly well 
suited to the present task have been developed in |2"4] . 

The outline of the paper is as follows: after introducing the multi-quark 
Lagrangian in Section 2.1, the stability conditions of the vacuum are dis- 
cussed in Section 2.2. In Sections 2.3 - 2.5 we use bosonization and heat ker- 
nel methods to transform the multi-quark into a mesonic Lagrangian, and 
extract the relevant contributions to the mass spectra and gap-equations. 
Section 3 is dedicated to the analysis of the pseudoscalar observables and 
includes a detailed discussion of decay constants with particular emphasis 
on the relation of our one angle approach to the two mixing angle analysis. 
Explicit formulae for masses and mixing angles are obtained in Section 3.4. 
In Section 4 are presented the characteristics of scalars, numerical results 
are given in Section 5 and conclusions in Section 6. Two appendices contain 
respectively the detailed derivation of the uniqueness of the solutions of the 
stationary phase equations and of the solution of the equations which yield 
the matricial coefficients relevant for meson mass terms. 
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2. The model 



2.1 The multi-quark Lagrangian 

We discuss the system of light quarks u, d, s (Nf = 3) with multi-fermion 
interactions described by the Lagrangian 



Quark fields q have colour (N c = 3) and flavour indices which are suppressed. 
We suppose that four-, six-, and eight-quark interactions £ 4(? , C 6q , £ 8q are 
effectively local. Likewise, they are constructed from local quark bilinears, 
like the scalar S a = q^ a q, or the pseudoscalar P a = qi^^X a q "currents" . Such 
bilinears have the appropriate quantum numbers to describe mesons. This 
approximation corresponds to the task considered: we want to obtain, after 
bosonization, the tree level effective meson Lagrangian, with local vertices 
and local meson fields, and relate the coupling constants and masses of such 
a Lagrangian with the parameters of the quark model. Meson physics in the 
large N c limit is described by a local Lagrangian of this type [25] . 

The global chiral SU(3)l x SU(3)r symmetry of the Lagrangian (JTJ) at 
m = is spontaneously broken to the £77(3) group, showing the dynamical 
instability of the fully symmetric solutions of the theory. In addition, the 
current quark mass m, being a diagonal matrix in flavour space with elements 
diag(m u , rrid, m s ), explicitly breaks this symmetry down, retaining only the 
reduced SU (2)j x U(l)y symmetries of isospin and hypercharge conservation, 
if m u = m d ^ m s . 

The leading order (in N c counting) Lagrangian of light mesons and the 
corresponding underlying quark Lagrangian must inherit the U(3)l x U(3)r 
chiral symmetry of massless three-flavour QCD. In particular, it was argued 
that in the large N c limit of QCD with three massless quarks the pat- 
tern of spontaneous chiral symmetry breakdown is uniquely fixed: the chiral 
U(3)l x U(3)r group, under some highly plausible assumptions, necessarily 
breaks down to the diagonal U(3). In accordance with these expectations 
the short-range attractive U(3)l x U(3)r symmetric NJL-type interaction 



can be used to specify the corresponding part of the effective quark La- 
grangian in channels with quantum numbers J p = + ,0~ [21]. The ma- 
trices acting in flavour space, A a , a = 0, 1, . . . , 8, are normalized such that 



£ e ff = qiij^dft - m)q + C 4q + C 6q + C Sq + ■ ■ ■ ■ 



(1) 



A 9 = — [{qKq) 2 + {qiisKq) 2 



(2) 
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tr(A a A&) = 25 a b- Here A = y 1 1, and A&, k = 1,2, ... ,8 are the standard 
SU (3) Gell-Mann matrices. It is well-known that such four-quark interac- 
tions lead (for some values of the model parameters) to the formation of a 
quark condensate, which is invariant under the vector subgroup U(3) and 
thus breaks chiral invariance of the ground state in accordance with the re- 
quirements of three-flavour QCD. 

The 't Hooft determinantal interaction is described by the Lagrangian 

C 6q = «(det qP L q + det qP R q) (3) 

where the matrices Pl,r = (1 T 7s)/2 are chiral projectors and the deter- 
minant is over flavour indices. This interaction breaks explicitly the axial 
Ua(X) symmetry, lifting the degeneracy of r] and rj' meson masses (in the 
chiral limit), and violates Zweig's rule j2Hj due to flavour mixing. It affects 
also the scalar singlet and octet states pushing down the mass of the SU(3) 
singlet. 

The large N c behaviour of the model is reflected in the dimensionfull 
coupling constants, [G] = M~ 2 , [k] = M~ 5 , which count as G ~ 1/N C , 
k ~ 1/Nc f ■ As a result the NJL interaction (J2J) dominates over £ 6(J at large 
N c , as one would expect, because Zweig's rule is exact at N c = oo. These 
couplings have opposite signs: G > 0, k < 0. 

The eight-quark Lagrangian which describes the spin zero interactions 
contains two terms: £$q = C&) + Cfq [18 , where 

4? = %i UP R q m )(q m PLqi)f = § [tr(5 - iP)(S + iP)f 

= 9 i(Sl + Pt)\ (4) 
4? = 1 Q92[(qiPRq m )(q m PLqj)(qjPRqk)(qkPLqi)} 

= f^tr [(S - iP)(S + iP)(S - iP)(S + iP)\ 
lb 

92 



— tr (S 4 + P 4 + 4P 2 S 2 - 2PSPS 



16 

= — [dabedcde (S a SbS c Sd + P a PbP c Pd + 2S a SbP c Pd) 
o 

+ 4:facefbdeS a SbPcPd] ■ (5) 

Here the trace is taken over flavour indices i,j = 1,2,3; the matrices S, P 
are given by Sij = S a (X a )ij = 2q j q i , P„ = P a (X a )ij = 2^(^75)%. The 
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f a b c are the well-known totally antisymmetric structure constants: [A ffl , A&] = 
2if abc \ c . The d abc are totally symmetric quantities: {A a , A;,} = 2d abc \ c . C Sq 
is a U(3)l x U(3)r symmetric interaction with OZI- violating effects in C&j . 

The eight-quark interactions £$ q are the lowest order terms in number 
of quark fields which stabilize the vacuum state of the model. We restrict 
our consideration to theses terms, because in the long wavelength limit (or 
in the case when the multi-quark correlators create a hierarchy) the higher 
dimensional operators are suppressed. 

Since the coupling constants G,K,gi,g 2 are dimensionful, the model is 
not renormalizable. We use the cutoff A to make quark loops finite. The 
regularization procedure (Pauli - Villars) is standard and can be found, for 
instance, in our paper [21], where the regularization function is introduced 
to define the coincidence limit of the Schwinger - DeWitt representation for 
the real part of the quark-loop effective action. This method is used for our 
calculations of mass spectra in the following. 



2.2 Stability conditions for the vacuum 



The eight-quark forces stabilize the vacuum [T^]. To clarify the meaning 
of this statement, consider, for the sake of simplicity, the effective potential 
V(M) which one obtains as a result of bosonization of these multi-quark 
vertices in the chirally symmetric limit (m = 0) and in the one-quark-loop 
approximation 



V(M) 



h 2 / 27 
(l2G + Kh + —Xh 2 



16 

3N C 

167T 2 



/ M 2 
M 2 J (M 2 )+A 4 ln 1 + 



(6) 



with A being an ultraviolet cutoff in the quark one-loop diagrams, and 

' A 2 \ 



Jo(M 2 



A — M In 1 + 



M 2 



(7) 



The dependence on the variable h is defined by the stationary phase equation 



Q 

M + Gh + — h 2 + -Xh 3 
16 4 



A = gi + -02 



(8) 



as a function of the model parameters and the argument M. 
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We start the discussion of the effective potential with the standard case 
of four-quark interactions, where the curvature of the potential at the origin 
and the sign of the coupling G of the interaction fully determine the existence 
of a globally stable system, which can occur either in the Wigner-Weyl or in 
the phase of spontaneously broken chiral symmetry. Often the word "insta- 
bility" is used in connection with the transition from this symmetric to the 
spontaneously broken vacuum at a critical value of GA 2 . This is not what is 
meant when we say that the vacuum is unstable. As we hope will be clear 
after the discussion presented in the remaining of this section, the instability 
we refer to is an essential pathology of the vacuum, present in the model with 
combined four and six quark interactions: we show that it is crucial for the 
stability of the vacuum that the stationary phase equation JSJ) possesses only 
one single real root when higher order multi-quark interactions are present. 
We argue that the enlarged system with six-quark interactions fails in this re- 
spect and that eight-quark interactions are necessary to stabilize the vacuum. 
Fig. 1 will illustrate the various stages of the discussion. 

In a world without six and eight-quark interactions, K,gi,g 2 = 0, one 
obtains from Eq. (jHJ) that h = —M/G and, as a result, the potential V(M) 
has the form of a double well, if V"(0) < 0, see Fig. Id, i.e. if r, the following 
combination of model parameters, 

N C GA 2 

* = -b->i- <») 

This inequality expresses the fact that chiral symmetry is spontaneously bro- 
ken 5 , producing in the massless case, (m = 0), the degeneracy of a nonet of 
Goldstone bosons, and showing the presence of the C/a(1) problem. This vac- 
uum state is globally stable, because at large values of \M\ another inequality 
is fulfilled 

Of* 9 Jl/T2 

V(M) ~ — h 2 (M) = — > (M^±oo) (10) 

at G > 0. In this particular case the stability of the ground state (G > 0) is 
already guaranteed by Eq. Q. 

It is worth noting that higher order multi-quark interactions will not 
change condition (jHJ) as long as Eq. (jHJ) has only one real solution. The 
reason for this is very simple. If this equation has only one real solution, it 

5 The Wigner-Weyl phase appears for < r < 1, i.e. V"(0) > 0, see Fig. la. 
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Figure 1: The effective potential V in the SU(3) chiral limit, as function of the 
quark condensate related variable h. Upper and lower panels are classified by the 
value of r = Nc ^t , related to the curvature of the effective potential at the origin. 
Each panel shows the typical form of the potential when one adds successively to 
the four-quark (left), the six- (middle) and eight-quark interactions (right). The 
metastable cases (b) and (e) are obtained in the mean field approximation Vmf 
(the stationary phase approach leads to an unstable vacuum, without any local 
minimum, |17j). 

is valid to expect that h(M) = —M/G + 0(M 2 ). Since 2n-quark vertices 
contribute to V(M) as h n ~ M n , the value of ^"(0) is entirely determined 
by terms of the second power in M, i.e. by the four-quark interaction (n = 2) 
only. Therefore it is tempting to describe the general situation (when higher 
order multi-quark interactions are included) by the same inequality (jHJ), since 
it will dictate the behavior of the effective potential in the neighbourhood of 
the origin in full agreement with the leading order result; this is illustrated in 
Fig. 1, the upper panel for r < 1, the lower one for r > 1. For that one must 
find however a way to reduce the number of real roots of the corresponding 
stationary phase equation to one. 

What is wrong with several roots? Let us consider the system which 
includes four- and six-quark interactions G, k ^ 0, gi, g 2 = 0. In this case 
the quadratic Eq. (jHj) has two solutions, both being real for M > AG 2 / k. It 
follows then that the stationary phase method leads us to the gap equation 



9 



which contains the sum of these solutions 



+ h& + ^£ J (M 2 ) = 0. (11) 



7T 2 



The sum does not depend on M, because h« + /i( 2 ) = -16G/«, and this 
gap equation misses the trivial solution M = 0, corresponding to the chiral 
symmetric vacuum. One sees that a simple addition of the 't Hooft inter- 
action to the four-quark Lagrangian affects so violently the trivial solution 
and as a matter of fact the whole effective potential, which gets unstable 6 , 
that, apparently, we must get rid of the problem which appears as soon as 
the stationary phase equation has more than one real solution 7 . 

Since a quadratic equation never has only one real root, we are pushed to 
increase the order of the equation by including eight-quark interactions. One 
obtains in this way the cubic equation (jHJ). This equation has only one real 
solution h(M), which changes smoothly in the open interval — oo < h < oo 
being an isomorphic and monotonic function of M, when one restricts the 
choice of parameters to 

G>\(^)\ A>0. (12) 



A V24 

For this case the constituent quark mass M fulfills the gap equation 

h(M) + ^fjo(M 2 ) = 0, (13) 

related with the potential V(M) (see Eq. (JBJ)) which is bounded from below 
(see Figs, lc and If). 

Eqs. (fT2|l replace the previous requirement G > (see Eq. (fTOJl ). They 
must be fulfilled to guarantee the global stability of the system. The first 
inequality is new and plays for the enlarged system the role of Eq. fllOj) : in 
the case of four-quark interactions only, the linear stationary phase equation 



6 This point has been considered in detail in |17| . 

7 Here we would like to stress that Figs, lb and le, related with the addition of six-quark 
interactions, are obtained within the mean field approach, which leads to the effective 
potential © taken at A = and considered as a function of h ^j. We identify it with 
Vmf(1i). The dependence M(h) is given by Eq. (JHJ). This is a one-to-one mapping 
h —* M, where h ranges along the interval — oo < h < oo. The local maximum at positive 
h = —8G/k on both figures corresponds to the point where the regular (at K — > 0) solution 
changes to the singular . VmfQt) is unbounded from below, as — > oo. 
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had automatically only one real root, here the values of couplings must be 
fixed correspondingly to ensure the existence of only one real root. The 
second inequality is a direct analogue of G > 0. 

Eq. is still relevant to the case and is responsible for the behaviour of 
V(M) in the neighbourhood of zero, as mentioned before. In Figs, lc and If 
the stabilizing effect due to the addition of eight-quark interactions is shown. 
Note that they change radically the potential only at values of h > —8G/k, as 
compared to the cases lb and le, calculated in the mean field approximation 
(see also footnote 7), affecting little the other branch of the potential, where 
h ~ hS 1 '. In particular the value of h where the global minimum of both 
potentials occurs in the spontaneously broken phase, is negative. Since at 
the quark one-loop order h is proportional to the quark condensate |22|, one 
is inclined to believe that by fixing the model parameters through it, this 
will finally lead to similar numerical values for all observables which depend 
in a stringent way on the value of the condensate. For those observables, the 
calculatinos in the metastable mean-field approximation and in the globally 
stable case considered with inclusion of the eight-quark interactions will not 
differ much. 

Another interesting aspect of this simple analysis of the SU(3) limit of the 
effective potential is the possibility of existence of multiple vacua, illustrated 
in Fig. 2. For r < 1, i.e. in a region where the four-quark interactions 
alone lead to the symmetric Wigner-Weyl phase, the inclusion of the 't Hooft 
six-quark intearctions can induce spontaneous symmetry breaking for some 
critical values of the coupling parameter k. This new vacuum can coexist 
with the trivial vacuum. This has been discussed previously in [HUj, now 
we are able to confirm its existence within the general stability conditions 
imposed for the vacuum, given by the inequalities in (|12j) . In the presence of 
non- vanishing values for the current quark mass, the minimum at the origin 
shifts towards the physical region of negative h. The phenomenon of multiple 
vacua has been addressed in several other approaches to the description of 
the QCD vacuum [HI1 IH2 E3 • 

The most general case is obtained if the SU (3)l X SU (3)# chiral symmetry 
of the quark Lagrangian is broken down explicitly by the non-zero values of 
current quark masses m. Then the inequalities (fT2"j) must be replaced by the 
following ones |18j 




(14) 
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V 




Figure 2: A closer view on the effective potential V of Fig. lc. Depending on the 
strength of the six-quark coupling k, with remaining parameters G, A, A fixed, the 
symmetric Wigner-Weyl phase of the four-quark potential (Fig. la) may coexist 
with the spontaneously broken phase induced by the presence of the 't Hooft term. 
This double vacuum exists within the global stability conditions (|12)l. 



The last constraint can be used to make a large N c estimate for g\. Indeed, 
we know that G scales as 1/N C , k ~ 1/A^ and, therefore, conclude from the 
above inequality that g\ cannot scale as or smaller. On the other hand, 
this eight-quark interaction is an additional (to the 't Hooft determinant) 
source of OZI-violating effects and thus it cannot be stronger than the 't 
Hooft interactions, i.e., g\ ~ or ^ ess - These reasonings show that 

1/A^ < gi < 1/iVj. We would expect #2 to be in the same interval, although 
inequalities (fb"fj) are not appropriate to prove that. 



2.3 Semi-bo sonized Lagrangian 

The multi-quark Lagrangian can be presented in the bilinear form with 
respect to the quark fields. The details can be found in fTJ Ej- With- 
out using any approximations one obtains the following vacuum-to-vacuum 
amplitude of the theory 

Z = / VqVqY[Va a Y[V(f) a exp (i / d 4 x£ q (q, q, a, 4>) ) 
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-oo 

X 



+00 

J Y[^s a Y[T>p a exp^z J d 4 x£ r (cr, 0, A; s,p)J , (15) 



-00 



£ 3 = g(«7 M ^ - M - a - 2 75 0)g, (16) 
C-r = s a (a a + A a ) + p a (f) a + — (s„ + Pa 



where 



2 

+ ^ ^Uc^a (SfeS c - 3p b p c ) + y( S 1+ Pi 

+ — [dabedcde {s a S b S c S d + 2s a S b p c p d + PaPbPcPd) 
o 

+ 4:facefbdeS a SbPcPd] • (17) 

The bosonic fields o = a a \ a and (p = (p a X a are the composite scalar and 
pseudoscalar nonets which will be identified later with the corresponding 
physical states. The auxiliary fields s a and p a must be integrated out from 
the effective mesonic Lagrangian C r . The quarks obtain their constituent 
masses M = M a \ a = diag(M M , M d , M s ) due to dynamical chiral symmetry 
breaking in the physical vacuum state, A a = M a —m a . The totally symmetric 
constants A a b c are related to the flavour determinant, and equal to 

1 

3! 

Some useful relations for A abc can be found in jSHj. 



A-abc — —\ e ijk e mnl\^a) im{^b) jn{^c) kl- (18) 



2-4 Stationary phase approximation for Z 

The functional integrals over auxiliary variables s a ,p a in Eq. (|15|) can be 
calculated approximately within the stationary phase method. For that one 
should first find all real stationary phase trajectories s a l = s a (cr,4>), p a l = 
p a {c, <p) given by the system of equations 

^ = 0, |^ = 0. (19) 

OS a dp a 

We seek these solutions in form of expansions in the external mesonic 
fields a a , <p a 

s% = h a + h a l b ] (7 b + h a l b \a b a c + h a 2 l<f) b <f) c + ... 

Pa=h a 2 Ub + h^ b a c + ... (20) 
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The coefficients h$ depend on the coupling constants G, k, gi, g 2 and quark 
masses A . The higher index coefficients h$ are recurrently expressed in 
terms of the lower ones. The one-index coefficients h a are the solutions of 
the following system of cubic equations 

A a + Gh a + —A abc h b h c + -^h a h\ + -^d abe d cde h b h c h d = 0. (21) 

Thus, the problem is reduced to a finite set of algebraic equations which span 
the nonet space of U(3). These equations may be considered as an example 
of response equations, previously studied in [33j for the purpose of gaining 
some understanding about the possible origins of strong SU(3) breaking. The 
response h a is fully determined by the couplings G,n,gi,g2 and the mean 
field A a , which plays the role of an external force. In accordance with the 
pattern of explicit symmetry breaking the mean field has only three non-zero 
components, with indices a = 0, 3, 8. 

In order to solve Eqs. (j21j) . one has to find whether there exists an 
intersection of a number of hypersurfaces. The important question of com- 
pleteness of the system (J21j) is examined in Appendix A. This yields h a = 
for a = 1,2, 4, 5, 6, 7. Thus the system reduces to three coupled equations 
to determine ho, h^, h$- This task has been solved in [TH]. At this stage, as 
we have already discussed, one has to find conditions (see Eq. (JHj)) which 
ensure the one-to-one mapping A a <-> h a . 

The next two equations following from (j!9)l determine the coefficients 
h®, h a 2 } in (JHJ 

I (g + y hfj 5 ar + ^A abr h b + g\h a h r 

+ ^ (2d abe d rce + dcbedrae) h c h b \ h^J = -8 as , (22) 



+ (Zfabefrce + d cbe d rae ) h c h b ^ h^J = —5 as . (23) 

Corresponding solutions are given in Appendix B. 

This procedure can be easily extended. Equating to zero the factor at 
any independent field combination in l)19|) . one obtains an equation which 
determines one of the coefficients in fl2()l) . 
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On the other hand, these equations are useful if one wants to find the 
projection of the Lagrangian C r on the stationary phase trajectory (|2(Jj). For 
that one should rewrite them in a more convenient form, using that 



' h h a h 2 b - y ( 

-jh a h 2 b + 3y< 

+ 92 face frdeKhdhr ■ (25) 



(h^)-Jh r = Gh a + 2A a - ^h a h 2 b - ^d abe d rce h b h c h r , (24) 
-{h i2) )arh r = 3Gh a + 2A a + 3^-hahl + 3^-d abe d cde h b h c h 



In particular, solutions of Eqs. (|2T]) - (j2"3|) define the first three coupling 
constants of such Lagrangian, i.e., one can show that 

C r -> C st = h a a a + l - h { 2a a a b + ~ h^ a fe + 0(field 3 ). (26) 

Since the system of equations (fT^j) can be solved, we are able to obtain 
the semi-classical asymptotics of the functional integral over s a and p a in Z. 
If parameters G,n,gi,g2 belong to a range where the system has a unique 
real solution, the calculations are straightforward. In particular, one has 
the following result which is valid at lowest order of the stationary phase 
approximation 



J YlDsaYiVpa exp (i J d 4 xC r (a, 0, A; s,p) 

— oo a a 

~ exp (i J d 4 x£ st (a,0)) (h -> 0). (27) 



2.5 Integrating quark fields in Z 

To obtain the effective Lagrangian of the model we should integrate out quark 
fields from Eq. ()15)) . This is a well studied part of the calculations and we 
restrict ourselves to several general remarks here. 

The one-quark-loop effective action can be computed in euclidean space- 
time, the chiral invariant part of the result, W^cr, (/>], is given by the modulus 
of the quark determinant 

W q [<T,<p] = ln\detD E \, (28) 

where De stands for the Dirac operator in euclidean spacetime, namely 
De = ilfx^ii — M — a — z75</>. The quark determinant is a complicated 
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nonlocal functional which can be approximated in the low-energy regime by 
the Schwinger - DeWitt asymptotic expansion (221 122] • The presence of a 
noncommutative (with respect to the bosonic fields a and 0) mass matrix M 
requires a more delicate treatment of this term in comparison with the stan- 
dard approach, where M is supposed to commutative with the fields. The 
corresponding technique has been recently developed |24< and applied to the 
case considered here in j2H]- We refer to these papers for necessary details 
(see Section 3 in j2H]), although we present the result, because we need it in 
the following. 

The heat kernel expansion used is 

W,M] = -|^ (29) 

where coefficients 6j for the case with isospin symmetry are 

bi = -Y, 6 2 = ^ + ^|A 8 y, ... (30) 

Our following result is based on these two terms of the series. It is the lowest 
order approximation, because the usual kinetic term of the collective fields 
is contained in 6 2 , and we truncate the series exactly after this term. The 
part of & 2 with A us = M 2 — M 2 is absent from the standard Seeley - DeWitt 
coefficient a 2 . This is one of the new features of the approach, which follows 
from the noncommutativity of the constituent quark mass matrix M. 

The trace in Eq. (j2Uj) should be taken over colour, flavour and four-spinors 
indices. In Eq. (j3*Uj) Y is used for 

Y = *7„(<V + i l5 d^4>) + a 2 + [M, a] + 4> 2 + i^[a + M, </>]. (31) 

The factors Jj are given by the average 

1 r 



h = - 

3 



2J i (M^) + J i (M^) . (32) 



and represent one-quark-loop integrals. In the considered approximation we 
need only to know J (M 2 ) (see Eq. (J7|)) and 

J ^M 1 + w)-vTW- < 33) 
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We are using here the proper time regularization scheme. 

Thus, the integration over quark fields yields the second part of the ef- 
fective Lagrangian (the first part is given by Eq. ([26))) 



f~-q * f~"hk f-'tad ^kin ^~"m ^int- 



(34) 



The tadpole term, C ta d, is 



C 



tad 



127T 2 



[M u (3/ - A us h)(a u + a d ) + M s (3/ + 2A us I 1 )a s 



(35) 



The kinetic term, Ckin, after continuation to the Minkowski spacetime, 
requires a redefinition of meson fields to obtain the standard factor in front, 
i.e., 



C 



kin 



NJi 

167T 2 



tr 



(<V) 2 + (d, 



1 



-tr 



where 



_(2 n ~.Q> 



«k, 9 2 



4tt 2 

The contribution to the mass Lagrangian is given by 



(36) 
(37) 



Cr 



NJt 



47T 5 



Nh 



127T 2 

+2(2M 2 + M s > 2 + (M 2 + 5M> 2 + (7M 2 - M 2 )er 2 
+ (M U + M S ){M U + 2M> 2 . + (M s - M U ){2M S - M u )0 2 } 



{A us [2V2(3(7o(7 8 



(38) 



where we assume that the indices i and / range over the subsets i = 1,2,3 
and / = 4, 5, 6, 7 of the set a = 0, 1, . . . , 8. Thus we have 



a 



2n 7r + (n 



°^ 2 4> 2 f = 2(K + K-+K K°), 
i° )\ a 2 f = 2(K* + K*-+K*°KZ°) 



(39) 



3. Pseudoscalars: masses, mixings and all that 

3.1 Symmetry and currents 

One can simply obtain the conserved (or partially conserved) currents of the 
local theory by using the variational method of Gell-Mann and Levy |3*o] . 
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For instance, the infinitesimal local chiral transformations of the quark fields 
in C e ff (see Eq. (JTJ) are 

Sq = i( a + 7 5 /5)g, 5q = -iq(a - 7 5 /3), (40) 

where the small parameters a = a a 4f and (3 = /3 a 2 are Hermitian flavour 
matrices. Then, according to the Gell-Mann - Levy formula, one obtains the 
standard vector V® and axial- vector A ^ nonet quark currents 



ya 



5C 



eff 



S(da a ) 
and their divergences 



_ A a 
<77 M y9, 



AO, 

A 1 



5£ 



_ A a 



0M£ 



5£ 



<5a< a 



- q [m, X a ] q, 
l -qi5{m,\a]q 



+ z'5 o0 v6k (det qP L q - det qP R q) 



(41) 



(42) 



(43) 



Transformations (|4()jl induce the correlated change in the flavour space of 
collective fields 



5a 



>R : 



a, afi + Mg 



-1 



z a, 



+ {PAr}, 

+ Mg- 1 }. 



(44) 



The quark Lagrangian, C e ff, is approximated by the effective bosonized La- 
grangian C bos 

£eff £-bos = C-st + C-hk + • • • , (45) 

where dots correspond to all omitted terms due to the approximations made. 
Therefore, one can obtain the currents again, but now they will be written 
in terms of meson fields. Indeed, one has 

! [({o-fl + Mg~\ d^n) - {d^ R , <M) A a ] + 0(b 3 ), (46) 



A a 



4 



-tr 



K A* 



I 

"4 tr 



a R + Mg- L ,d^a R + [<f> R , d^n] A a + 



(47) 



These currents obviously depend on the order where the heat kernel series is 
truncated, because the coefficient b 2 and higher ones contain derivatives. The 
symbol 0(bs) shows that currents (|46J) and (}4~Tj) have been obtained from the 
Lagrangian Chk based on two terms of the asymptotic series ([29)1 . namely b\ 
and b 2 ■ 
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3.2 Decay constants of pseudoscalars 

Let us calculate matrix elements of axial-vector currents 

(0|^(0M(p))=*7°V (48) 
Using Eq. (J46j) one derives 

y-oo _ 2M M + M s u _ ^, 22 _ ^,33 _ M u 

3g 9 

P 44 _ ^55 _ ^66 _ r77 _ M u + M s 



2g 



f ss M u + 2M S ^ f08 _ f80 _^ M u -M r 
3# 3^ 

It is not difficult to relate these abstract values with the experimentally 
measured decay constants of physical pseudoscalar states P(x). For instance, 
the weak decay constants of the pion (f n ) and kaon (f K ) are defined by the 
corresponding isotopic components of the axial current Aj^ t2 and -4^ +l5 , i.e., 

(0\A^ i2 (0)Hp)) = iyftUpv (0\A 4 + i5 (0)\K(p)) = iV2f KPfl , (50) 

and, therefore, one finds 8 

, M u M s + M u 

Jtt = , IK = ~ • (51) 

9 2g 

These are the Goldberger - Treiman relations at the quark level. 
The decay constants in the rj — 7/ system are defined 

(0\AI(0)\P<P)) = ifpPr, (« = 0,8) (52) 
where P = 77, 77'. Each of the two mesons has both, singlet and octet compo- 

IlGIllS 

U° R \ _ ( cos9 p -sm6 p \( V '\ ( . 

{(/>%) { sm6 p cos6 p )\ V J- {06> 

This orthogonal rotation diagonalizes the kinetic and mass terms in the me- 
son effective Lagrangian. The mixing angle 8 P will be calculated in Section 
3.4. Consequently, from Eq. ()52|) one obtains the 2x2 matrix {fp} 

SfaX - ( f" fv\-( COS ^ ~ Sin ^ ^ I ) M 

Upl ~i4 gj_~v^gp cos ^ A/ 80 f 00 )' { } 



3 We use the normalizations /„ = 92.42 ±0.26 MeV, and f K = 113.00 ± 1.03 MeV jSHj. 
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This matrix depends on four independent parameters g, M u , M s , 9 p . An alter- 
native parametrization has been considered in pT7| I38j. where two constants 
/o, /s and two angles i?q> $8 specify the matrix fp. 

{fp}=( f I™t* ~/° Sin f° V (55) 
UPJ ^ / 8 smtf 8 /ocos^o / v ; 

There is a straightforward correspondence between the parametrization 
of Kaiser and Leutwyler and the model predictions. Indeed, one finds 



(fs) 2 = if*) 2 + (#) 2 = ^ (M 2 + 2M S 2 ) , (56) 

(/o) 2 = (/J) 2 + (#) 2 = ^ ( 2M * + M *) • (57) 

The formulae for the relations between mixing angles are given in the end of 
this section. 

As one would expect, the model predictions agree well with the general 
requirements of chiral symmetry following from chiral perturbation theory 
(ChPT), although the results differ already at lowest order. For instance, we 

3 3g* 

One can show that the second term on the r.h.s. is of order (m s — m u ) 2 and 
therefore must be omitted at lowest order of ChPT. The rest of this formula 
is a well known low-energy relation which is valid in standard ChPT. 

The mixing angles i^^o are small and $8 7^ #o- We have for their 
difference 

/ 8 /o sin (ft - O ) = + = ~ (M s 2 - M 2 



(fs) 2 = + VVls :T U> ■ (58) 



(59) 



2^ ^ 2 ^ V 2 " , lf w ^ 
6# 



Thus, $8 7^ $0 due to the SU(3) flavour symmetry breaking effect. Again, 
this result agrees with the ChPT formula, if one notes that the last term is 
a higher order contribution. 

The singlet decay constant f is equal to 
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In the //-extended version of ChPT there is the OZI-rule violating term in 
the effective Lagrangian, which contributes as f%Ai to the r.h.s. of Eq. (|6U|). 
We have instead the term (M s /M u — l) 2 /6 ~ Ai. Of course, in our case 
the origin of this contribution is related with the SU(3) flavour symmetry 
breaking. 

Let us now express the mixing angles i?s, in terms of one mixing angle 
6 p and quark masses. One has 

tan tf 8 = -± = tan \9 P - arctan — - — ] , ( (> I ) 



tan^ = - A. = tan ^0 p + arctan -^jr 1 ^ ) • (62) 



It follows then that 



tf 8 = V -arctan (v 7 ^) , (63) 
i?o = ^ - arctan , (64) 

where ip = 9 P + arctan a/2. Similar formulae have been obtained in jSH] by 
Feldmann, Kroll, and Stech. One should not confuse the angle 9 P , defined 
by the rotation ()53|) and contributing to the matrix fp as it is shown in Eq. 
(|54p. with the naive identification $ 8 = $ — @p discussed in the literature in 
connection with the one mixing angle problem. Our consideration is perfectly 
consistent with the two mixing angles approach. 

3.3 Strange - nonstrange basis for decay couplings 

The axial-vector currents can be taken in a different basis, namely, we shall 
consider now the nonstrange A^f and strange A s ^ currents 

The singlet <p R and octet <f) 8 R fields are also rotated to the new basis 

(66) 
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The corresponding matrix elements are easily calculated 

(0|^(0)|^(p)) = ip^ , <O|^(O)|0»> = ip^ . (67) 

The physical states P = 77, 7/ are the mixtures of the nonstrange and 
strange components, this follows from Eqs. and (jbbp 



cosip sintp \ (V 
— sin ip cos ip J \r)'. 



(68) 



where the angle ip — &p + arctan y/2 ~ 9 P + 54.74°. 

Next one can find the couplings describing decays of physical states in 
the hadron vacuum 

(0\^(0)\P(p))=if P p„ (i = ns,s). (69) 

The result can be represented in a way which is similar to the one of Leutwyler 
- Kaiser J40| 

ff i X _ ( f™ fv \ _ ( fnsCostf ns -/ s sintf s \ 

UP) [fl? flf) \fns^ns / S COS^ s ) ^ 

Our calculations show that 

rns M u , rs M s 

/ ns = COS ip, f* = sm V, 

9 9 

f ns=Ml sini)j f s i = M± cosij . (71) 

9 1 9 

It follows that the basic parameters f ns , f s , t? ns , $ s of the matrix {/}>}, being 
expressed in terms of model parameters (in the approximation considered), 
are 

fns = — = A, fs = —, Ip = $ns = $8- (72) 

There is a direct relation between a common mixing angle fi ns = fi s and the 
OZI-rule which has been discussed in p>9~]-|4*T]. 
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3.4 Mass formulae and the mixing angle 9 P 

The gap equations are an essential ingredient to obtain mass formulae of 
pseudoscalars. One obtains them equating to zero the tadpole contributions 
from the Lagrangian Cbos, see Eq. (j4*51) 



iV 

h u + -\M U (3/ - A us h) = 0, 
N 

h s + -\M S (3/ + 2A us h) = 0. 



(73) 



The stationary phase equations ()143j) taken in the isospin limit {m u = rrid) 
have been also used to obtain the mass spectrum 9 



' g 2 m u 



m 



GM U l + uo s + p + T u 



9 ( m u + m s 



K 



1 



G \M U + MJ 1 + u u + p + t uu + t s , 

2 



T,. 



m 



9_ 
2 



(A + B =p y/(A - Bf + AD 2 



(74) 
(75) 
(76) 



Here 7/_ =77,77+ = 77'. 

There is a mixing between the (0, 8) states in the multiplet. The sym- 
metric mass matrix M p is given by 



Mr 



9_ 
2 



>>Wk) 



A D 
D B 



where we have 



, „ _ K K 2(1 + p) -uj s + t uu + t s 
+ Mi, M, GdetiV(2) 



A-B 
D 



1 

3 



h s i 8u u + U s + T ss - T u 
GdetNM 



+ 



+ 



3 \M U M, 



Gdet M 2 ) 



(77) 

(78) 
(79) 

(80) 



The non-diagonal term of the matrix vanishes in the SU (3) flavour symmetric 
case (m u = = m s ), otherwise D ^ 0. This matrix is diagonalized by an 



3 See Eq. Ijl48(l for our notations of u>i, p and n- 
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orthogonal transformation (|53|) to the states (rj, rj') with masses given by Eqs. 
(I76j) . The singlet-octet mixing angle is 

ta n26 p = — -. (81) 

It is easily seen from Eq. (|8*T|) that the mixing angle 8 P is equal to its 
ideal value: tan(2# id ) = 2y/2, i.e., 6 id ~ 35.26° and ^ = 90°, if « = 0. As a 
result we have 10 rj oc </> ns and rf oc s , thus the eight-quark interactions have 
no influence on the flavour content of rj, rf without the 't Hooft term. 

Consider now the rj — rf masses (J75J) presented as follows 



m\ T =m 2 K + Q 1 T VKr - ml - Q 2 f + 2Q§ . (82) 
The independent functions Qi, Q 2 , Q3 are equal to 

g 2 (2uJ u + UJ s + T uu + T ss -2T us ) 2 

Ql ~ 2GdetiV(2) (M s -M u ) 

.2 



29 ( u « — t us )(uj s — (jj u + r uu + t ss — 2t us ) 
G(l + p + uj u + t uu + t ss - t us ) det AT( 2 ) 

Q2 ~ 2GdetiV(2) + ^~ M ^ 



(83) 



_l g 2 (q> s - qju + t us - r ss ) 

C(l + p + uj s + t uu ){\ + p + u) u + t uu + r ss - r us ) ' 

Qs = GdetlvW ' (85) 

In the large-iV c approximation {Q\,Qz,Qs — > Qi°, Q\° 1 Qz°) these func- 
tions are related to each other and to the ghost coupling A^, introduced by 
Veneziano in j^j, 

2 -Q[° = 2Q L 2 ° = Qs° = ^f = §- c (large N e ), (86) 
where ui ~ 1/N C is a leading order contribution of u u , or u s . We have for X 2 



2 = _ kNc_ /M\ 



(87) 

N c ^oo 

10 Note that the orthogonal transformation l|53() is written for k ^ 0, what corresponds 
to A — B > 0. If k = 0, one has the opposite inequality A — B < 0, and, as a consequence, 
one should replace in Eqs. Ij53(l and (|68|l the fields 77 •«-> J]'. 



24 



With these specific values of Q's our expressions for the masses of rj, 7/ mesons 
coincide with Eq. (34) of Veneziano work in jH]. Moreover, the Witten - 
Veneziano formula for the mass of 7/, 

ml / +ml-2m 2 K = X 2 v = —^x(0)\ Y M, (88) 

J 7T 

which is obtained in the large A" c limit of QCD for non-vanishing quark 
masses, relates the 7/ mass with the topological susceptibility in pure Yang 
- Mills theory x(®)\ym- Consequently, Eq. (|H7j) yields 

x(0)|yM = ^(^) 3 (large AQ. (89) 

The interesting feature here is the explicit demonstration that the eight-quark 
forces contribute to Eqs. (|83J ) - (j8l)j) in such a way that the only dominant term 
in (|88|) is still the 't Hooft interaction, even though the Pij, Tij may formally 
be of the same l/iV c -order as u>i. 

One can use the most recent lattice calculation of the topological sus- 
ceptibility in 02]: x(0)\ym = -(1-33 ± 0.14) x i(T 3 GeV 4 to find Q\ 
0.156 GeV 2 . Additionally, we have at leading A^-order 



tan », = 2V2 ml _ m , - ^ -+ 2V2 ml _ m , _ gQ , , (90) 

obtaining approximately a mixing angle 9 P ~ —14°, and masses m, q ~ 
516 MeV, m r] r ~ 1077 MeV. These numbers reflect the general picture pre- 
sented in Section 5 rather well. 



4. Scalars: masses and the mixing angle 9 S 

4-1 Mass spectrum of the scalar nonet 

The masses of the scalar nonet: a (/ = 1), Kq (I = 1/2), (I = 0), are 



2g 2 (u s - r u 

"" ' T ' G[(l + p + 2T uu y-(uj s -T uu )i} 



< =< + ±K + TWn ■ "V. -V, 7 _ , (91) 



m k = ml + 4M„M S + ^ - - - _ ^ , (92) 

m ^ = |- ^ + S zp ^ - £) 2 + 4D 2 ) . (93) 
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To obtain the last formula (|93|) the mass matrix of isospin singlets <j° r and 
a 8 R in the Lagrangian (JHS|) 



a*. = £m,«S)(£ v B ){°h, (94) 



where 



i , R _ K^.h^ NJi ( M 2 , »a\ , 2 + ^ + 41 + 3(7-^ + ^) 



h u 

— + 


h s 

— + 


K 


h s 


M u 


M s 


( K 


h s 




M s 



3GdetN^ 

/ u u 

V = V2 





3GdetNW y 
has been diagonalized by an orthogonal transformation 

(95) 

/ \ — fill's woi's / \U ftl 

with the angle given by 

TD 

tan20 s = — — -. (96) 

Let us apply formula ()96|) to the case in which the 't Hooft determinant 
is neglected. We find that the mixing angle 8 S ^ 9\& at k — 0. The reason is 
that the terms which are proportional to p us in T> and A — B have different 
coefficients. Therefore the /J" meson has an admixture of the strange com- 
ponent and correspondingly the f$ meson has an admixture of nonstrange 
quarks due to the eight-quark interactions 11 with coupling g\. Such admix- 
tures explicitly violate the OZI rule in these scalar channels. 

We must notice that the singlet-octet splitting of scalars is more sensitive 
to the eight-quark interactions, as opposed to the pseudoscalar case con- 
sidered above. We can gain some understanding of this by writing slightly 
less explicit formulae for scalars. For that we turn again to the large N c 
arguments, postponing the exact calculations till Section 5. 



1 We have / oc (uu + dd) and Jq oc ss at k = g\ = 0. 
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4-2 The 1/N C consideration 

The 1/N C expansion is usually a good approximation for hadrons. If we 
accept this idea, we can deduce from the above formulae a clear qualitative 
picture of the role played by the eight-quark forces in the mass spectra of 
scalars. Our starting point are the following expressions 



< = K + AM'i + ^{uj s -r uu ) + ... , (97) 



2q 2 

m 2 K * = m\ + \.\I„.\L + 4t K -r us ) + ... , (98) 



2 - m 2 K + (M s + M u ) 2 + Q 1 



T V (A^ + Q 2 ) 2 + 8(A K7T + Q 3 y , (99) 

where A Kn = m 2 K — m 2 , and ellipses are used to denote omitted terms of 
l/N 2 order and higher. The functions <2i, Q2, Q3 are 

Q 2 

Ql= 2G ^ ~ Us ~ 2 P~ ( Tuu + Tss + 2Tus ^ + ■ ■ • ' ( 10 °) 
Q 2 

Q2 = -Q2 + ^ (4w„ + 2cj s + 2p uu - p ss + 8p us ) + . . . , (101) 



Qs = -Q2 + ^ ( 2w. - y + 2 Puu - p ss - p us j +... . (102) 

The contributions from the different r's are exactly canceled in Eqs. (J101J) - 
()102j) at this order. Moreover, note that the SU(3) breaking corrections 
have a higher order in 1/N C . For instance, uj u ,uj s ~ l/iV c , but the difference 
uj u — uj s ~ ^/N 2 . The same is true for other flavour dependent functions in 
Eqs. (jl00|) - (jl02|) . Thus, one has to lowest order in l/N c 

\LO n,LO nrpLO nt?LO s^LO n ( n,LO , rpLO\ r\LO 



Q^ = Q% u -3Et; u -2EZ u , Q^ = 3(Q^ + Et u ), Q^ u = 0, (103) 



where the eight-quark contributions E\ and E^° , namely 

,2, 



rpLO 9 Puu 



E. 



LO _ 9 T uu 



2 ~ G 

large N c 



(104) 



large N c 



are proportional to the strengths ~ g\ and g 2 correspondingly. 

Next, let us try to understand in simple terms the hierarchy inside the 
nonet. It is easy to see from Eqs. (j9T|) - (j9l?|) that m r < m ao < m K * <m f +. 
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This is in an agreement with the result of analysis |43J . Let us do some crude 
numerical estimates. We have at leading order 12 



m 2 K * - m)j = Q L 2 ° ± y (AjfTT + Qk°) 2 + 8A^, 
2 - A K7T + 4M U (M S -M U ) } 



ni 



m„ 



A Kn = 2M s (M s - M u ). 



(105) 
(106) 
(107) 



First, we conclude that eight-quark forces are probably unimportant for 
the difference m 2 K * — m 2 . 

Second, we use the ratio fx/ fir = 1-22 to find M s /M u . Indeed, it follows 
from Eqs. (jBTjl . (fHfi jl . Ij57j) . (fHS jl and (ffi4 ]l that 



If M s 
- 1 + — 

2 V Af„ 



A 

U~\A Ml, 



$o — ^8 — arctan 



fn \A Ml 





(M s 


M u \ 




\M U 


MJ_ 



(108) 
(109) 
(110) 



One easily finds M s /M u = 1.44. Numerically, this yields / 8 — 1-31/tt, fo = 
1.17 U, ^0-^8 = 19.5°. 

Third, we use phenomenological values m n ~ 138 MeV (averaged over the 
isotopic triplet n°, tt^) and nix — 495.7 MeV (averaged value for isotopic 
duplet K+, K°) to obtain A Kn = 0.227 GeV 2 , and m\» - m 2 o ~ 2.39 A Kn = 

0.542 GeV 2 . 

In order to make more progress, we need a further dynamical input. For 
that we identify the quark- ant iquark ao-state of the model with the known 
I G (J PC ) = 1~(0 ++ ) resonance ao(980). This resonance is often considered as 
KK molecular-like bound state jH]. The four-quark nature of the ao(980) 
meson is also widely discussed in the literature (see, e.g., the recent paper [13] 
and references therein), where ao(980) is a compact KK state. The extended 



12 The last Eq. (|107|l follows from our general result Eqs. (|7^jl - l(75|l in the pseudoscalar 
sector for the difference 



2M S (M S - M u 



0(1 



T"as ) 
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molecule case does not exclude that the core part of the wave function may 
be dominantly qq (SI- The four-quark picture is however based essentially 
on the MIT-bag model, thus representing an alternative approach to the 
problem. 

Now, by using that m ao ( 980 ) ~ 980 MeV, one derives rriK* — 1226 MeV. 
We suppose that this state may be identified with the wide I(J P ) = |(0 + ) 
resonance AT* (800): m K * {800) = 797 ± 19 ±43 MeV, T = 410 ±43 ±87 MeV, 
reported in [4*7] : see also |H]-jS2]; where further support for this low lying 
state is given. 

Before we calculate the masses of the two states from Eq. ()105j) . con- 
sider the singlet-octet mixing angle 9 S which can be written very compactly 
in the large N c world, viz., 



If the 't Hooft interaction term Q\° dominates over E\° in (this is 
possible when g x ~ l/iV c 5 ; in this case E\° ~ l/iV c 2 < Q\° ~ 1/N C ), one 
easily finds that 9 S ~ 21° (we have used here the estimate Q\° = 0.156 GeV 2 
obtained before Eq. (jHOJ)). 

One immediately derives from Eq. ()105j) the masses of the singlet-octet 
mixed states frf: m f - ~ 300 MeV and m f + ~ 1407MeV. 

There is a strong cancellation in the formula for the /o~-mass 



m). ~ m 2 K . - 3Q%° - \J (A K7T + 3Q%°y + 8A^ 

~ (1.502 - 0.468 - 0.946 = 0.088) GeV 2 . (112) 

As a consequence the result is very sensitive to the parameters of the model. 
We shall see later that the best fit of the pseudoscalar channel leads to the 
value rrij- ~ 550 — 750 MeV. Therefore, the lowest mass scalar meson, J'q, 
is identified with / (600). The Particle Data Group assigns to this resonance 
the mass m /o(600 ) = 400 - 1200 MeV, and the width T = 600 - 1000 MeV. 

The state Jq agrees with the state /o(1370): mf ( 1370 ) = 1200 — 1500 MeV 
(r = 200 — 500 MeV). The lower state /o(980) with the same quantum num- 
bers is too far away from the value following from Eq. (J105)) . thus our esti- 
mate shows that /o(980) may be not a member of the scalar quark- ant iquark 
nonet considered 13 



13 



It is notorious that the low-lying scalars are still the subject of many studies. Our 
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It is known that the U(1)a breaking produces an opposite in sign mass- 
squared splitting between the octet and the singlet for the scalar and pseu- 
doscalar mesons. This can be embodied in the approximate sum rule |l3*l 154] 

ml + m 2 - 2mi + ml+ + m 2 - 2mi* ~ 0. (113) 

'' '' Jo Jo o v ' 

The r.h.s. results from the exact cancellation between the 1/N C order terms 
induced by the 't Hooft interaction (see, for example, Eq. ijBB]) ). One can 
see that eight-quark interactions may contribute to the sum rule already at 
l/N c order, if the coupling gi counts as Indeed, one obtains 

m% + m 2 v - 2m 2 K + mj+ + m). - 2m 2 K * = -QE^° + O J . (1 14) 

This term has a negative sign, decreasing the sum m 2 + m 2 + . Let us note 

A) Jo 

that in this case g\ lowers the value of m^- and due to the fine tuning effect 
in Eq. ()112j) the octet-singlet splitting grows with increasing gi in the scalar 
nonet. Thus the above sum rule is a good illustration of the possible impact 
of the eight-quark OZI violating forces on the scalar mesons. 



5. Numerical results 

We collect the results of the exact numerical calculations for the mass-spectra, 
mixing angles and quark condensates in three tables. Input is denoted by a *. 
Table 1 contains the seven parameters of the model m u , m s , A, G, k, g±, g^. 

Sets (a, b, c) and (d, e, f), are each a block for which we compare the effect 
of the new parameters gi, g2 as follows. In the first line of each grouping we 
set #1, gi to zero, and fit four of the remaining parameters (m u , m s , G, k) by 
fixing m^, rax, fn, fx] in the first set A is fixed through /q~, in the second 
through r] f . The reason to fix the empirically not well known mass of /J" is 
that it is the most sensitive to changes of the parameter gi (we remind our 
discussion of Eq. (|112|) ). By fixing it, one reverts the situation and is able to 
detect the effects of g%, gi on the other observables. In the last set one sees 
that by fixing the mass of 77', it is Jq that monopolizes the value of g\. 

conclusion agrees with some other results 021 > but we are aware that the point 
requires an additional analysis of scalar decays to be definitive. In different approaches, 
based on a coupled channel analysis, the pole position may also be affected by closed 
channels, see e.g. 
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Table 1 

Parameters of the model: m u , m s (MeV), G (GeV~ 2 ), A (MeV), re (GeV -5 ), g x , g 2 
(GeV~ 8 ). We also show the corresponding values of constituent quark masses M u 
and M s (MeV). 





m u 


m s 


M u 


M s 


A 


G 


— K 


9i 


92 


a 


5.2 


161 


302 


486 


934 


7.18 


1122 


0* 


0* 


b 


5.5 


175 


325 


523 


896 


8.78 


774 


1000* 


0* 


c 


5.4 


173 


322 


519 


900 


8.57 


822 


1000* 


-132* 


d 


6.1 


189 


372 


646 


839 


12.16 


1082 


0* 


0* 


c 


6.1 


189 


372 


646 


839 


11.28 


1083 


1500* 


327.24 


f 


6.1 


189 


372 


646 


839 


8.92 


1083 


6000* 


327.24 



For all sets, except of course the case g\ = g% = 0, the stability conditions 
dHJ) are fulfilled. 

Tables 2 and 3 contain the results to be compared with the phenomeno- 
logical data. The first set shows drastic effects of the parameter g x on mass 
spectra: going from (a) to (b) the rj' mass is reduced by 40%, getting close to 
its empirical value m v t = 957.78 ± 0.14 MeV, while the rj mass gets smaller 
by 7%. The overall effect is a reduction of the gap between these two states, 
which can be translated into the smaller value of the parameter re in (b), as 
compared to the one in (a). 

The rj — rf splitting can be illustrated by the formula (|82|) which takes 
the form 

2 2 8(m? K - ml) 2 + 3c q 

™>r, = m 7T{ 2 2^ ' 115 

' 9(777y - TOq) 

where = ^(4m 2 K — m 2 ) is the Gell-Mann - Okubo result for the r^-mass. 
The remainder originates in the repulsion of r\ and rj' and represents an SU(3) 
breaking effect of second order. The coefficient c q depends on the Q's given 
by Eqs. (JHl-dHHI), namely, c q = 2(m 2 K - m 2 7r )(Q 1 - 3Q 2 ) + 3(Q| - Q\ + 2Q 2 ). 
Formula ()115j) for c q ^ extends the Veneziano result [Hj (see Eqs. (34) 
there) by including the SU(3) breaking corrections stemming from the 't 
Hooft and eight-quark interactions. 

Numerically m = 565 MeV is just a little bit larger, to be compared 
with the phenomenological value = 547.30 ± 0.12 MeV. The Witten- 
Veneziano correction (the second term ~ (to^ — to 2 ) 2 ) is related to the 
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topological susceptibility 14 and is about four times larger than it is required. 
Considering set (b), we obtain ~ 496 MeV. This value is now corrected 
by the six and eight-quark contributions collected in c q . Unfortunately they 
work in the same direction and do not improve the low result for rj. Finally 
one finds the value = 486 MeV. This is the general tendency in all sets. 

Table 2 

The masses, weak decay constants of light pseudoscalars (in MeV), the singlet- 
octet mixing angle 9 p (in degrees), and the quark condensates (uu), (ss) expressed 
as usual by positive combinations in MeV. 





m 


7T 


m K 




Try 


U 


h 




i 

— (uu) 3 


-(ss)^ 


a 


131 


I* 


494* 


525 


1761 


92* 


120* 


-1. 


246 


210 


b 


131 


S* 


494* 


486 


968 


92* 


120* 


-12 


242 


199 


c 


m 




494* 


493 


1023 


92* 


120* 


-10 


242 


200 


d 


m 




494* 


476 


958* 


92* 


116* 


-14.4 


233 


184 


c 


m 


** 


494* 


476 


958* 


92* 


116* 


-14.4 


233 


184 


f 


m 


V* 


494* 


476 


958* 


92* 


116* 


-14.4 


233 


184 



Other effects are seen in the reduction of all scalar masses, except the 
input one, by 25% for a , 14% for Kq, 9% for f£. In set (c) the effect of 
g2, keeping the same g\ as in set (b), is seen to increase all masses again. 
It has in this case a negative value, which is also allowed by the stability 
conditions. The condensates and scalar mixing angle remain comparable in 
(a, b,c), the pseudoscalar angle 9 P increases in absolute value from set (a) 
to (b). Its value quoted in Table 2 is around 6 P = —12° being in agreement 
with estimates of Veneziano jH] and more recent calculations in [55j. 

This angle is directly related by Eqs. (|BHJ)-(|EIJ) with two mixing angles 
seen in the singlet and the octet components of the decay constants. We 
find, for example, for set (c): $$ = —21.6°, $o = 3.5°. 

There is some ambiguity about the definition of quark condensates if the 
chiral symmetry is explicitly broken by bare quark masses mj. The values 
given in Table 2 are obtained by the subtraction of the expectation value of 
qtqi in the perturbative vacuum from its expectation in the true vacuum: 

(Mi) = 2 ( h iU^o - h *lA l= o) > ( 116 ) 
14 One can find details, in particular, in the first reference of 
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which is the definition used in [TU] . Let us recall the recent update of the light- 
quark condensate at a scale of 1 GeV: (qq)(l GeV) = —(242 ± 15 MeV) 3 , 
where qq = (uu + dd)/2 represents the isospin average of the non-strage 
quarks The flavour breaking ratio is known to be (ss) / (qq) = 0.8 ± 0.3 

inm 

The second set (d, e,f): we fix the rj' mass to its empirical value, and 
chose also fx closer to experiment, keeping m n ,mK, f n as in all other cases. 
In (e) we obtain g% through the mass of ao(980) and take g\ arbitrarily. In the 
pseudoscalar sector this lowers slightly the value of the //-mass and changes 
a bit the mixing angle 6 P . 

Table 3 

The masses of the scalar nonet (in MeV), and the corresponding singlet-octet 
mixing angle 9 S (in degrees). 





m a (980) 


^X*(800) 


m /o(600) 


m /o(1370) 


S 


a 


1262 


1347 


600* 


1436 


16 


b 


945 


1150 


600* 


1309 


15 


c 


980 


1176 


600* 


1326 


21 


d 


993 


1217 


754 


1391 


25 


c 


980* 


1204 


691 


1374 


23 


f 


980* 


1204 


559 


1362 


23 



The main effect is visible in the /o(600) mass: a further increase in gi, 
set (f), decreases further its mass, leaving the remining observables almost 
unaffected. The repulsion between the two isosinglet levels caused by the 
eight-quark interaction thus lowers the value of m/ (6oo) by about 200 MeV. 

There is no marked effect on the -K"q(800) state, which continues to lie 
above the ao(980) mass. The large N c result ()106j) protects the inequality 
rriK* > TTi ao f° r the members of the quark- antiquark octet, which finally holds 
in the general case. 

To summarize, the effect of eight-quark interactions on the mass spec- 
trum, vacuum decay couplings, and mixing angles is relatively small as long 
as general properties of the QCD vacuum (the values of the quark conden- 
sates and the topological susceptibility) are correctly reproduced. 
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6. Conclusions 



The role played by eight-quark interactions in the long wavelength limit of 
QCD has been addressed in a systematic way. A full understanding of its 
impact on the vacuum and properties of the low-lying spin zero mesonic 
spectra has been achieved. As an important by-product also the results 
associated with the well known four and six quark Lagrangians due to Nambu 
- Jona-Lasinio and 't Hooft, on the body of which the eight-quark terms are 
attached, are classified and presented according to stability criteria of the 
effective potential (reviewed and illustrated in Section 2.2), the large N c 
counting scheme, approximate sum rules for meson masses, U^(l) and flavor 
SU(3) breaking terms, and OZI-rule violation. This "dissection" allows not 
only for a complete understanding of the full result, (i.e. in the leading order 
stationary phase approximation to the functional integral), for which we also 
give analytical expressions, but also to compare with other relevant works in 
the field. These discussions and respective formulae accompany every main 
derived step within the full result and serve as a "hitchhiker's guide" to the 
mesons in the multi-quark "galaxy". 15 

An issue of much interest is the two-angle analysis of the rj — rj' mixing 
and its relation with the standard one-angle diagonalization. We clarify 
some existing confusion in the literature by deriving in full detail in Sections 
3.2 and 3.3 the connection and equivalence of the two methods within the 
Lagrangian considered. 

The masses and splitting of the complete result for 77 — 77' system are 
conveniently cast in the form of Eq. (jll5|) . which separates in a transparent 
way the leading Gell-Mann - Okubo contribution, the Witten - Veneziano 
correction and the second-order SU(3) flavor breaking corrections due to six 
and eight-quark terms. These latter ones have a positive sign for the param- 
eter sets which yield good fits for the remaining pseudoscalar observables 
and therefore add a small correction to the already large Witten - Veneziano 
term, thus yielding a larger splitting than the empirically observed. 

Concerning the scalar sector we show that there exists a mass hierar- 
chy within the model considerations which is not conform with the present 
understanding of the empirical results. 

These drawbacks, being the result of an exhaustive and consistent study, 
clearly indicate that effects not considered, such as meson loops, higher orders 

15 Borrowed from "The ultimate hitchhiker's guide to the galaxy", by D. Adams. 
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in the heat kernel expansion, and confining forces, might be at work. 

We view therefore the main role of eight-quark forces considered as folows: 
(i) they are of vital importance for the stability of the ground state built 
from four and six-quark interactions. They restrict the choice of the cou- 
plings G, k, gx, g2 to the rather narrow window of combinations given by Eq. 
(I14|) . This is important, since combinations outside the allowed range can 
at instances even yield a better spectrum for the pseudoscalar mesons alone. 
It would be an erroneous result, attributing minor importance to the cor- 
rections of the kind not considered, (ii) They help in understanding the 
effects caused by the OZI- violating terms with coupling strength gx, which 
affect quite strongly the splitting of the f$, /q" scalars, mainly pushing down 
the lower state, due to the strong cancellations reported in Eq. (|112|) . (iii) 
They may be also of importance in decays and scattering, not considered so 
far. (iv) They give a clear indication that an hierarchy of multi-quark in- 
teractions, with dominance of lower ones, is present. This corroborates with 
recent lattice calculations jSZI of gluon correlators, where the lower ones also 
dominate. 
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Appendix A. Response equation (12 ip 

An algebraic system of nonlinear equations (}2"T]) has, in general, more than 
one admissible solution. Our specific problem requires only the knowledge 
of all isolated real roots. In order to find them suppose that the system has 
at least one real solution, i.e., the set {h a }. Let us now manipulate with 
the equations at hand in such a way that we can guess some values from 
this set. If our guess is correct, the system reduces to a smaller one, which 
finally can be solved. It may happen however that the smaller system is 
incomplete and therefore has a continuum of solutions. Such cases have to 
be excluded, because they are possible if and only if some of the components 
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A a are constrained (see below). It should be recalled that A a is a set of 
independent variables (for different values of a) which will be fixed only later 
with the help of the gap equations. Thus any restriction on A a at this stage 
leads to an internal contradiction and the corresponding solutions must be 
rejected. There is another reason to exclude a continuum of solutions, namely 
the stationary phase method cannot be applied for such a case. 

After these general remarks, let's consider important details. Eq. (12 lj) 
yields the following set 



phi + z/(/i 4 /i 6 + h 5 h 7 ) 

ph 2 + v(h 5 h 6 - h 4 h 7 ) 

A 3 + ph 3 + V - (hi + h 2 5 -h 2 6 

ph 4 + v{h x h % - h 2 h 7 ) 
ph 5 + v{h x h 7 + h 2 h 6 ) 
ah 6 + v{h x h 4 + h 2 h 5 ) 
ah 7 + v{hih^ — h 2 hi) 



0. 
0. 



h 2 7 



0. 

o, 
o, 

0. 



0. 



A 8 + i(p + cx) h + ^ 



hl + h 2 2 + h 2 3 -^ (h 2 4 + h\ + hl + h 2 7 



(117) 

(118) 

(119) 

(120) 
(121) 
(122) 
(123) 

0, (124) 



Ao + St* + |^/f (3h 2 - h 2 a ) + ^=h s [h\ + h 2 2 + hl-fj 



hi + hi) (vsh 3 - h 8 ) - (hi + hi\ (V3h 3 + h 



92 

2V2 



+W - [hi (h 4 h 6 + h 5 h 7 ) + h 2 (h 5 h 6 - h A h 7 )} = 0, 



(125) 



where /1, u, p, a and £ are defined as follows 
H = G + ^_ (^ -h ) + ^ (gi + g 2 ) h\ + |/i (ho + 2y/2h 8 ) , (126) 



" = - ! u A l^h , £ = G+l(g 1 + 2g 2 ) h 2 a - ^h 2 , 



K 

16 



1 
2 



(127) 
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p = G- 



a = G- 



K 

16 



K 

16 



-ho — /la H 1^ 

3 3 ^3 



j +^ (<7i + <7 2 ) hl+^ho (h + v^a - V2h 



SJ, 

(128) 



y^|/io + h 3 + +^ + <? 2 ) /i 2 + |/i (/iq - VQh 3 - v/2/ig) • 



Due to these definitions one has 

p - a = v (V3h 8 + hfj , p - p = v (\/3h 8 



fi- 



(130) 



We start from the observation: if p ^ 0, one can multiply Eqs. (jl2()|) - 



(131) 
(132) 
(133) 
(134) 



by p and, using Eqs. (|117|)-(jll8j). obtain 




hi 


pp - v 2 (hi + hfj 


= 0, 


h 5 


fip - v 2 (hi + hfj 


= 0, 


h 6 


pa - v 2 (h 2 4 + hfj 


' =0 ' 


h 7 


pa — v 2 (hi + hfj 


' =0. 



This leads to four possible alternatives: 

1. /i ^ 0, hi = h 5 = h$ = h 7 = 0. 

2. p j£ 0, h 4 = h 5 = 0, a = 0. 

3. p 0, h 6 = h 7 = 0, p = 0. 

4. p ^ 0, /ip = i/ 2 + hfj , pa = z/ 2 + hfj . 

If p = 0, then Eqs. (|117|) - (jll8|) now imply other four alternatives: 



5. 


p = 


o. 


v = 0. 








6. 


V = 


0. 




/ii 


= h 


= 0. 


7. 


V = 


0. 




/i 6 


= h 7 


= 0. 


8. 


p = 


o. 


i/^0, 


/; i 


= h 5 


= h 6 = 



h 7 = 0. 



(135) 
(136) 
(137) 
(138) 



(139) 
(140) 
(141) 
(142) 



Cases 2 and 3, as well as cases 6 and 7, are correlated. The existence 
of such correlation becomes clear if one notes the invariants of Eqs. (|117|) - 
(|125|) under the substitutions /i 4 <-> h 7 , h 5 h 6 , h 3 <-> —h 3 , A 3 <-> — A 3 . 
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In particular these substitutions change p <-> a, leaving p, and v without 
changes. 

Case 1. In this case Eqs. flll7|) - ()118j) yield h\ — h 2 — and one obtains a 
set of three equations to determine ho,ti3,hg. To make further progress one 
should switch to the flavour basis 0, 3, 8 — > u, d, s, where one has 

' Gh u + A u + ^ h d h s + ^ h u {h 2 u + h 2 d + h 2 s ) + ^ h 3 u = 0, 

lb 4 I 

Gh d + A d + ^ + j + h 2 d + h 2 s ) + | h\ = 0, ( 143 ) 

GTi, + A s + ^ + ^ ^ s (^ 4-^ + ^)4-^^ = o. 
lb 4 Z 



This system has been studied in our work 

Case 2. Again Eqs. ([TT7il-(fTT8l yield h x = h 2 = 0. Then Eqs. and 
fT^l give A 3 = VSA 8 (with the use of Eqs. l|T3n|l and (7 = 0). Since the 
variables A a are supposed to be independent, one concludes that this result 
leads to an apparent contradiction. It may be also noted that the pair of 
quantities (A 3 , A 8 ) define in general a plane in the nonet space, if neither 
quantity vanishes. The found correlation means that only one axis is defined, 
in which case the system is not complete and the solution is underdetermined. 

Case 3. As before Eqs. flHTjl - fTTHj) yield hx = h 2 = 0. Then Eqs. (fTT3j) 
and ()124[) with the use of Eqs. (JI30J) and p = give A 3 = -v^Ag. One 
comes anew to the abovmentioned contradiction. 

Case 4- From Eqs. ()138j) one obtains p(p — a) — v 2 {h\ + h 2 — h\ — h\). 
Next, note that one can use Eqs. ()130|) to rewrite the result as follows 



2ph ?J + v (hl + hi- hi- h 



0. 



Since v ^ this yields A 3 = 0. (The equality v = would lead here 
to p = a = 0, and finally due to Eqs. ()130j) to p = 0. The obtained 
contradiction proves that v ^ 0.) Eqs. ()117|) - ()118p . as well as equations 
fll38j) . imply that pa = u 2 {h\ + h 2 2 ). Using this result and Eqs. (jl3()j) . one 
can show that A§ = 0. It is clear that in this case the solutions suffer from 
the same kind of defects just mentioned above. 

Case 5. Since p, — v — 0, Eq. ()130J) gives p = a = 0. This yields A 3 = 
and A 8 = 0. This system is incomplete. 

Case 6. Eqs. (IT22J) -(IT2?1 yield a = 0. Since p, = a = 0, Eq. ([HDl gives 
^-3 + y/3h$ = 0. If he and h 7 are nonzero, then from Eqs. ()12()|) and ()121|) . 
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h\ + h\ = — > hi = h 2 = 0, since we only allow real solutions. As a result, 
one obtains 

A 3 - v^Ag = uh 3 ( h 3 + V3h 8 ) = 0. 



Case 7. Eqs. ffT20jl and (jT7f]l yield p = 0. Next, due to a = 0, Eqs. 
(I122j) and ()123j) reduce to + h 2 h^ = 0, /ii/is — h 2 h^ = 0, it follows then 
hi = h 2 = 0. Eqs. ()13()j) give /i 3 = v^3^8> that finally leads to the correlation 
A 3 = -v^Ag. 

Case 8. Eq. (ITPIJ) reduces to A 3 = 0. Eq. (IT2H) has the form 
A 8 + + h 2 2 + h 2 3 - 3hl) = 0. 

This case gives a new class of solutions. Nevertheless it can also be thrown 
out, as long as A 3 = 0, by the same arguments as just given. 

One concludes that only the case 1 leads to the isolated real solutions 
of Eqs. (PTTj) which must be taken into account. As has been shown in [TS] 
one can choose the parameters of the model in such a way that only one real 
solution appears. 

We would like to note that our investigation here recalls in many aspects 
the old results of Pais [34J, who studied the response equations with octet 
driving forces. He dealt with the octet space of SU(3). Eqs. (J2T]) have some 
more complicated structure, but the general conclusions remain true. Even 
the formal covariance property of the system under the transformation 



h = qo, h = -q 6 , h 2 = q 7 , h 3 = - (q 3 - V3q 8 ) , h± = g 4 , h 



2 



15 = 95, 



h = qi, h 7 = -q 2 , h 8 = -\ (qs + v 7 ^) , (144) 



in the nonet space is fulfilled. This transformation brings Eq. (|21j) in the 
form 

b a + Gq a + — A abc q b q c + —q a ql + -^d ahe d cde q b q c q d = 0, (145) 

where 

h = \ (A 3 - v^Ag) , 6 8 = -i(A 8 + v / 3A 3 ), (146) 
and all other b a = 0. 
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Appendix B. Solving equations (1221) and (|23|) in the isotopic limit 



As it has been shown in Appendix A, the system (|21J1 is actually reduced 
to three coupled equations (|143|) for three independent variables h u ,hd,h s . 
These variables are defined as follows h a X a = diag(/i u , ha, h s ), and easily 
related with ho, h% and h$, namely 

h u = —= (V2ho + V3h 3 + hi 
V3 



hd = — = (V2h - V3h 3 + ^ 



V3 

h s = J^(h -V2h s ). (147) 

If we choose the current quark masses suitably, m u = ^ m s , the 
flavour symmetry of the action is broken down to the isospin group SU(2). In 
this partial case one finds from Eqs. (|2"T|) that h a = for a = 1, 2, 3, 4, 5, 6, 7, 
i.e., only two components ho and h$ are generally nonzero. Thus, in virtue 
of Eq. ()147|) one concludes that h u = h d ^ h s . We shall suppose that h u and 
h s are known 16 . 

Let us solve Eqs. (|2"2"j) - (}2"3"j) to find the couplings h a ]j and h^ for the case 
with isospin symmetry. To represent the result we shall use dimensionless 
quantities 

gihihj g2hihj 
<*=16G' ^ = r - = ^' (148) 

induced by six- and eight-quark interactions with strength couplings k, g\ 
and g 2 correspondingly. We use also that p = 2p uu + p ss . 
The result for /iv* is 



h {1) 

a ab 



-s, 



ab 



G(l - lo s + p + 3t uu ) 

Sab 



G(l -u u + p + r uu + r ss + r us ) 
-JV# 



GdetNW 



(a, b = 1,2, 3), 

(a, b = 4, 5, 6, 7), (149) 

(a, b = 0,8). 



16 See ^H] for details, where a more general case, h u ^ hd ^ h s , has been considered. 
The case studied here is a straightforward consequence of that result. 
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where the 2x2 symmetric matrix A^ 1 ) has elements 



N, 



(i) 



oo 

(i) 



1 + - (Uj a ~ 4:UJ U + 10p uu + 7p ss - 8 p us ) + T uu + 2T s , 



a/2 

N m = N m = -g" k - w s + 2(p ss - 2p uu + p us ) + 3(t ss - T uu )\ , 



(1) 



1 + 3 ( 4 ^« + 2uJ s + 14p„„ + 5p ss + 8p us ) + 2t uu + r s 



and its determinant is equal to 

det iV« = 1 + u s - 2u 2 u + 4p + 3 [(1 + p){r uu + t ss ) + p 2 
+ u s (p ss - 2p uu + r ss ) + 3r 2 s + 6r us p^ 



For h^j one obtains 



/i 



(2) 
ab 



-6, 



ab 



G(l + U S + p + T uu ) 
-Kb 



(7(1 + ^„ + p + r UM + 



-N, 



(2) 



06 



I GdetA^ 2 ) 

where the 2x2 symmetric matrix A^ 2 ) has elements 



(a,b = 


1,2,3), 


(a, b = 


4,5,6,7), 


(a, 6 = 


0,8). 



ZV (2) 
iV oo 


1 

= 1+ 3 


(4cu u 


1 

- w s ) + p + - 


(j~uu + 2t ss ) , 


iv 08 


= Aff = 


_ V2 
3 


(u s — UJ U + T ss 


Tuu) , 


iv 88 


2 

= 1 "3 


(2u u 


1 

+ U s ) + p + - 


(2t uu 4~ T ss ) , 



and 



(150) 



(151) 



(153) 



det = (1 + p)(l + p + t uu + r ss - Us) - 2u\ + r ss (r uu - u s ). (154) 
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